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Abstract 

We study the one component model for four different lattice 
actions in the Gaussian limit and for the Ising model in the broken 
phase. Emphasis is put on the euclidean invariance properties of the 
boson propagator. A measure of the violation of rotational symme- 
try serves as a tool to compare the regularization dependence of the 
triviality bound. 
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1 Motivation and introduction 



There is ample evidence that for the lattice regularized $^-model the physical 
quartic coupling decreases logarithmically with increasing cutoff and that the 
continuum limit describes non-interacting boson fields. Thus the scalar sector 
of the standard model 

where the Higgs fields live in 0(4), has to be considered as an effective 
low energy model describing the leading cutoff effects and may be studied 
independently. 

The cutoff then represents new, possibly more fundamental physics which 
would replace the standard-model at higher energies. However, its influence 
is effective on all scales and manifests itself in the cutoff dependence of all ob- 
servables. In describing the experiments through a lattice model, the higher 
in energy we go the stronger we expect such cutoff effects to be. More and 
more terms would have to be included in the effective action in order to 
explain experiments with sufficient accuracy. Within the scaling region one 
expects to find this dependencies to be only weak, but this knowledge is 
mainly based on rough estimates and there are only very few explicit studies 
i I, ill- 

In this article we study the cutoff dependence with emphasis on the eu- 
clidean invariance properties of the boson propagator. As a simplified toy 
model we use the one-component <I>^ model, where one does not have the 
additional complication due to Goldstone modes. We expect that most of 
our qualitative results carry over to 0{N). 

One particularly interesting observable exhibiting a cutoff dependence 
as discussed above is the upper bound on the effective coupling A/jQ. This 
bound emerges at infinite bare coupling (the Ising limit of the model) as a 
consequence of triviality, and is related to the Higgs mass via the ratio 

R = \/W3 = niH/f. (2) 

In first order perturbation-theory and for the 0(4) model the W mass may 
be related to the vacuum condensate f^^'^'^ = {^r) ~246 GeV, and therefore 

"'^The subscript R denotes renormalized quantities throughout the article. 
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(0) provides an upper bound for the physical Higgs mass 0. Because this 
quantity is essentially nonperturbative, we need nonperturbative methods to 
study its cutoff dependence. 

In our further discussion we refer to different lattice regularizations only. 
In lattice formulations the cut-off dependence enters through the discretiza- 
tion of the action and the geometry of the lattice mesh. In particular we 
expect regularization dependent lattice artifacts for large values of the di- 
mensionless Higgs mass (low values of the correlation length). Neuberger 
et al.0, have studied these effects in the 0(A^)-model by systematically 
adding operators of increasing order in 1/A; they also studied the contribu- 
tion of similar higher order terms in a 0(4) model on an F4-lattice 0. 

Different lattice geometries correspond to different lattice regularizations 
and have been studied mainly in the context of the model with 0(4) symme- 
try, in connection with the triviality bound on the Higgs mass. First results 
in were followed by a particular careful study of the model on the F4 
lattice 1^. During completion of this article we learned of a detailed study 
of the Symanzik action [^]. 

Here we present a Monte Carlo study of the one-component $4-model in 
the broken phase. For the standard hypercubic lattice this model has been 
extensively studied both, analytically |jlO|] and in various Monte Carlo simu- 
lations |Tl|. We investigate it on altogether four different lattice geometries, 
introduced in the context of the 0(4)-model and report on a simultaneous 
analysis. We compare the rotational invariance (RI) properties of the boson 
propagator in the broken phase of the corresponding Ising model with those 
of the related Gaussian model. We use the violation of this invariance as a 
measure in order to estimate the regularization dependence of the triviality 
bound. 



2 Various lattice regularizations 
2.1 Lattice actions 

Consider the action for the one-component $^-model defined on a lattice. 
We study four different kinds of lattice-geometries and corresponding ac- 
tions, which we call C, C, -F4 and S. All four lattices may be embedded 
into the regular hypercubic one C if" of volume L^, with different cou- 
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plings to nearest, next-to-nearest and on-axis distance two neighbour spins. 
For convenience we present all results on this hypercubic mesh with lattice 
constant a = 1. The general action S may then be written in the form 



N 



(3) 



By x + ej, we denote the neighbour site in direction e^, which may be a linear 
combination of the cartesian basisvectors in 4 dimensions {e^,/i = 1, ...,4}. 

Table 1 summarizes our notation, in particular the number N of inter- 
action neighbours, the weight factors rjy and a normalization factor a that 
has been defined in such a way that the trivial continuum limit of the kinetic 
term is universal (with a leading term in the inverse momentum space 
propagator) . 

We denote the 'lattice field' by $ to distinguish it from the field $ in 
'continuum' normalization. The naive continuum action is reconstructed in 
the limit a ^ 0, L ^ oo with the substitutions 



mo 



l-2g 1^ 

l^Vv 

an a ^ 

[any 

C denotes the usual simple hypercubic lattice action with interactions 
between nearest neighbours along the directions of the cartesian axes. In 
comparison to C the C-action has additional couplings to the next- 
to-nearest neighbours in diagonal directions so that there are altogether 32 
interaction neighbours. -F4 is the nearest neighbour action on an F4-lattice 
[y, 0. It's embedding into a hypercubic lattice may be imagined by removing 
all odd sites of the grid. The symmetries for this case forbid terms like X^/^ 
in the inverse momentum space propagator and guarantee Lorentz-invariance 
to a higher order in the momentum cut-off (0(A~^)) than the first two hy- 
percubic lattice schemes (0(A^^)). The S denotes the Symanzik improved 
action which includes antiferromagnetic couplings to the the next-to-nearest 
neighbours along the 4 euclidean axes-directions. These additional terms are 
incorporated in order to remove the undesired 0(A~^) contributions, like for 
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i<4- The remaining corrections of 0(A ^) in the lattice dispersion relation 
guarantee tree- level improvement (cf. sect. 4). 



With 



the kinetic parts can be written 

4 



kin 



^kin 



'-'kin 



-2k ^x+e,^a 
,1=1 

4 

,1=1 

4 ,1-1 

^=1 1^=1 

-2k J2 i^x+e 

,1=1 V 



u=l 



x+Ef^+e^ + ^x+efj,-e„ j '^x ; 
^x+2eu 
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(5) 

(6) 
(7) 
(8) 
(9) 



Note that C^*^ does not couple odd with even sites. 

Quantization amounts to summing over all field configurations and taking 
expectation values with regard to the Gibbs measure exp {—S) where 
the partition function Z is the normalization factor such that (1) = 1. For the 
Gaussian version of the model this can be done explicitly, since the functional 
integrals factorize in momentum space. For the general theory one has to 
rely on expansion techniques or Monte Carlo integration. 



2.2 How to compare two regular izat ion schemes 

Comparing results from different actions obtained at finite cut-off is not 
straightforward, though. Different actions have different scaling violation 
properties. Choosing couplings such that the observed correlations lengths 
agree is not necessarily the best definition. Consider e.g. an effective action 
for block variables (which usually is more complicated). Comparing the 
results of the original and the block action at a point in their respective 
coupling constant spaces where the dimensionless correlation lengths agree, 
the block action will definitely have better continuum properties. It is "closer 
to the continuum limit" . 
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So one has to define what one means by the intrinsic ratio between the 
scale parameters in two different regularization schemes A and B. In an 
asymptotically free theory that issue can be resolved by adjusting the be- 
haviour of the /3-function in the weak coupling region |jT3|, |T0|, |1^. In our 
case in principle we can rely on renormalized perturbation theory. In fact, 
the two-loop expansion gives 

m^hysas = C,(AA^)-^=^/^?e-i/^i^«{l + 0{\r)} (10) 

where the subscript s denotes the regularization scheme (lattice action or 
geometry). Asymptotically, at small enough small A/j, ag/ag' = Cg/Cg' can 
be determined by lattice perturbation theory (cf. the results of llOl for the 
hypercubic lattice model). 

For a comparison of different effective theories at finite cutoff and not too 
small Xr we can no longer rely on that approach. It would be appropriate to 
choose another physical observable O, that may serve as a means to measure 
the distance from the continuum limit independently. Comparing results due 
to different actions Sg at the same values of O then provides an estimate for 
the regularization dependence. This observable could be for example a tttt- 
scattering cross-section, the a-decay width or the violation of euclidean RI 
of some quantities 0, |^, ^. Different operators O will usually have different 
cut-off dependence and their choice needs further motivation: There is no 
"best" observable. Quantities like the ratio R should be considered versus 
such a "distance" O to the continuum limit. Of course we could still use the 
inverse correlation length ruphygag, but in order to compare different schemes 
we have to rescale at each point such that O agrees. 

In this study for four different kinds of lattice actions we compute the 
envelope of R in the Ising limit {g = oo). In order to relate the respective 
results we use for O the RI violations of the inverse correlation length and 
the inverse 1-particle momentum space propagator. 

3 Rotational invariance (RI) 

The process of regularization by a lattice implies a violation of rotational 
symmetry for non- vanishing lattice constant. Only in the continuum limit 
do we hope to restore full euclidean symmetry. Thus the amount of violation 
of RI provides a measure of how close one is to the continuum. 
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In order to estimate the violation of Rl one can in principle refer to 
scattering cross sections. However, measuring these in computer simula- 
tions is notoriously difficult. For that reason one relied on tree level lattice 
perturbative calculations in order to estimate the leading violations for e.g. 
Goldstone-Goldstone scattering in the 0(4) model [|l^, |l^, H, |, H] where the 
Born term is proportional to the massive scalar propagator. In the one com- 
ponent model this would require the study of scattering of massive Higgs 
particles. 

In our approach we want check the violation of RI in the Higgs propagator 
That quantity is more accessible in the simulation. For simplicity we 



consider only 2-d slices of our hypercubic lattice. This has the additional 
advantage that the Monte Carlo data results from operators summed over 
the other two directions provide better statistics. 

The problem then is to quantify the RI of some function F defined on a 
planar grid. Specifying F to be a real positive function in the plane, with 
the symmetry property F(r, ip) = F(r, —ip) = F{r,(p+ |7r), fcGZ, ipE [— vr, vr) . 



Then F has a series-expansion of the form 

oo 

F{r, ip) = Y^ a„(r) cos{4rnp) 

n=0 



(11) 



For fixed r we define p>Q to be the directional angle for which F is minimal. 
For the inverse propagator our models always have ipo = 0; for the inverse 
correlation length p}Q = for C, C and = 7r/4 for i^4, S (in a d=2 slice). 
A measure for the deviation from the RI of F(r, p>) at the distance r can be 
defined by 



\ 



1 

2^ 



F{r,ipo) 



dp> 



\ 



;i2) 



n=l 



where 



A„(r) 



F{r,ipo) 



The last relation follows from the orthogonality of the cosine functions. The 
leading contribution to J\f may be estimated from the ratio ^p'^^^^ since in 

leading order J\f 



F(r,7r/4) 
F(r,0) 
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We can estimate the violation of RI either in real space or, by Fourier 
transformation, in momentum space. For infinite lattice size the momentum 
variables cover finite continuous intervals. In our discussion of the Gaussian 
model we utilize this property. 

In a straightforward application one could study the real space propagator 
for given fixed coupling k. However, it is clear that small and large distances 
(for periodic b.c.) trivially show strongest RI violation. In order to compare 
results at different values of the coupling (and the correlation length) we 
have to specify a typical value of | a; | or \p\ where we determine M for that 
coupling. We choose | p |= m{K,), a momentum value equal to the inverse 
correlation length in the direction of the cartesian axes. 

For all lattice actions considered we have calculated Af for 

i) the inverse momentum-space propagator G~^, denoted by Afc, and 

ii) the inverse directional correlation-length C,~^, denoted by A/"^ . 

In the continuum limit m — > we have A/" ^ 0. 

We determined A/" for the Gaussian model both analytically with Taylor- 
expansion and by use of conventional numerical methods. For the Monte 
Carlo results of our simulation of the Ising model we relied on suitable in- 
terpolation of the numerically determined momentum space propagator, as 
discussed later. 



i) Inverse momentum-space propagator 

We restrict the momentum-space to two dimensions by considering a particle 
which moves only in 2 of the 4 space-time dimensions with 4-momentum 
P = {Pt,Px,Py = 0,Pz = 0), with the notation pt = mcosy?, Px = msinyj, i.e. 
|p|= m. 

If the propagator is known explicitly, as it is the case for the Gaussian 
model, the simplest way to calculate Afcijn) is by the series-expansion of 
in the functions cos(4n<y9) (s denotes the lattice-action). The corresponding 
coefficients and A„ may then be expressed through Bessel functions and 
decrease rapidly. For n>l and p<l the ratio A„+i/An ~ lO"^". Therefore 
it is sufficient to consider only the first 2 coefficients [|16]. For propagator 



values given numerically one has to rely on an adequate interpolation or 
fitting procedure. 
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ii) Inverse directional correlation length 



The inverse directional correlation length ^^^{mslr) for the lattice- 
action s in direction r = (^i, f2, fs, f4), |f| = 1, is defined implicitly by 
the equation 

G;\K,\p = iC'r) = . (13) 

In that expression describes the asymptotic, exponential decay of the 
configuration-space propagator in direction r. In particular, along the main 
axes, e.g. ,^~"'^(m|(l, 0, 0, 0)) = m(K). Again we study the behaviour in a 2-d 
plane. 

Even in the Gaussian model the coefficients a„ and A„ cannot be rep- 
resented exactly, but one may express them through a power-series. The 
qualitative behaviour for m < 1 does not differ from the one of G~^. For 
numerically given propagators one may determine the directional correlation 
length from e.g. a fit along the specified direction (to the suitably inter- 
polated propagator). If the deviation from RI is correctly described by the 
equivalent Gaussian model, one may rely to a fit to the corresponding prop- 
agator. 



4 The Gaussian model 

Removing the quartic term in (^ we are left with the Gaussian model de- 
scribing a free scalar particle. The action is a positive-definite quadratic 
form, 

TV 

S = ^xQx,y^y ) Qx,y — ^x,y T]i^^x,y-ei, ) (l^) 

and it is possible to solve explicitly for quantities like the propagator. For 
the momentum-space propagator one finds 

N 

G{k\p)~^ = 2 1 ~ K'^r]^cos{p-eu) , (15) 

1^=1 

where the components of p assume values n = 0, . . . , L — 1 on an L"^ 
lattice with periodic boundary conditions. The configuration-space propaga- 
tor cannot be represented in similar closed analytical form. When 
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it is summed up over 3-space, however, it is possible to give an explicit ex- 
pression. For periodic boundary conditions x+Le.^ = x and if there are only 
ferromagnetic interactions (like in the cases C, C and the well-known 
result (cf. [|18|) is 



^($0$,) = -— ^ G U I = — — , ) e^P** 



L V ' L ' y 2(l-e-^'")sinhm 

(16) 

(Note that this is the propagator for fields in the continuum notation.) The 
relation between the lattice mass m and the hopping parameter k is defined 
implicitly by the pole of the momentum-space propagator in the complex 
energy-plane 

G-^ [k I im, O) = (17) 

with the solution 



smh— = -W . (18) 

2 2V a/t 

This implies the critical values for Kc = jj for these three actions (with 
corresponding values of N). 

For the case S, with the additional antiferromagnetic interaction term, 



(|T^ assumes the form 

8 - k(45 + 16coshm - cosh2m) = 0, (19) 

which has, for < k, < kappac = ^, two real solutions mi and 1712 obeying 
the relation coshmi + coshm2 = 8. Only mi approaches zero at the critical 
point. The other pole corresponds to a 'ghost' with a negative residue, ap- 
proaching 1712 — cosh~^(7) at the critical point. In the continuum- limit the 
corresponding physical mass goes to infinity and this state decouples. 
For this action ([T7| ) assumes the form 
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cosh 1712 — cosh mi 



_|_ (3-^i{L—t) g— '"12* _j_ g-m2(-L— <) 



1^1 — e"^'"!) sinhmi {1 — e~^"^'^) smh.m2 



(20) 



The closer to the phase transition we are, the more the physical mode domi- 
nates the propagator. Practically the ghost becomes noticeable only at small 
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distances t. E.g. for mi = 1 att = the ghost term contributes about 20% to 
the propagator, at t = 1 it is 4% , but at t = 8 its contribution has already 
decreased by several orders of magnitude and may be neglected. 

We now turn to the RI properties of the propagator. In the Gaussian 
model we have the advantage that we can work on both, finite and infinite 
volumes. We consider a two-dimensional submanifold, either in real space or 
in momentum space. 

Fig. la gives the inverse momentum space propagator G~^{p,(f) for the 
Gaussian models and a given value of m = 0.407, and in Fig. lb we plot the 
corresponding yj-dependence for p = | as an example. 

Following the above discussion and requiring | p | = m, in Fig. 2 we show 
lines of constant values of G~^{px = msmip,pt = m cos 99), as derived for 
infinite lattice size. For small m (small lattice spacing) rotational symmetry 
is restored. 

All this information may be summarized by plotting A/'g and Af^, shown 
in Fig. 3 for p < 1. In each case the RI is violated only up to a very 
small amount, expecially for F4, where A/'g,A/'^ < 0.0004 for that region 
of values. Even for the standard hypercubic action C where the largest 
violations appear, Af does not exceed the value 0.012. 

Generally in two or more dimensions is only calculable approximately, 
Taylor-series in the lattice- mass m^, 

C\rns\r) =ms[l + as{r)ml + 0{mt)] . (21) 

Both Rl-functions Mg and Af^ may be expanded 

Afpims) = A, {mly^ [1 + Bi^pml + 0(m^)] . (22) 

The leading coefficient is independent of F G {G~^, ^~^}, but it is not 
allowed to conclude that Ag is a universal constant 0. The same holds for 
js] it is a measure for the degree of euclidean symmetry restoration in the 
limit rris and has for C, C the value 1 and for F4, S the value 2. There 
are no O(a^) corrections for F4 and the tree- level improved Symanzik action 
S, which shows the higher degree of RI for these actions, one of the reasons 
for their construction. 
Postulating 

AfBiruB) =AfA{mA) (23) 
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we could rescale the corresponding definitions of tlie lattice constants and 
get asymptotic ratios 

^^v^ , ^ — ^ , (24) 
ac as 

for the Gaussian models. 



5 Results for the Ising limit 

The upper bound on the renormalized coupling and on the Higgs boson mass 
is realized for bare coupling g oo, the Ising limit. Thus it has become 
customary to study $^-models in that simple limit restricting the range of 
values of the field variable to ±1. 

In our Monte Carlo simulation we worked on lattices of size 16^ and peri- 
odic boundary conditions. Each of the 4 actions was considered at 4—6 values 
of the coupling parameter n > n^r in the broken phase, close to the critical 
point Kcr such that mG (0.2, 0.9). For the updating we used the cluster algo- 
rithm |T9[ (actions C, C and -F4) or the Metropolis method (action S, where 
we had problems with the cluster updating due to presence of additional an- 
tiferromagnetic interaction terms). In the determination of the propagators 
we utilized reduced variance operators ||20[. For each data point we consid- 
ered 40 000 (cluster algorithm) or 100 000 (Metropolis) update-measurement 
pairs, in addition to a suitable number of equilibrating updates. 

In this way we obtained the boson propagator summed over y, z directions 

G(/€|t, X) := ($0,0,0,0 E ^t,x,y,z) (25) 

y,z 

and its Fourier transform. We also determined the susceptibility 

^ = Y,G{K\t,x) = ^ . (26) 

In order to extract the required quantities = asUiphys and R, we used (0) 
and (^) for the Gaussian propagator 

G,(t) :=^G',(t,x) ^a(e-"^* + e-"^(^-*))+6 , (27) 

X 
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where b gives the square of the vacuum expectation value, ($)^. For our 
analysis we therefore use this value for | ($) |. The direct estimate of ($) is 
not possible due to tunneling between the two ground states. 

For the tree level improved action S we also add another term to ( ]27|) , 

^,^^-m'(m)t ^^-m'{m){L-t)^ ^ (28) 

where m'{m) is assumed to be given by the relation cosh m+cosh m' = 8 which 
is valid for the corresponding Gaussian model. 

The maximal relative statistical errors for the fit parameters amount to 
Am ^ ±2-3%, Aa ±1-2%, Ab ^ ±0.1% for results due to cluster 
updating and measuring. For the conventional updating and measurement 
(case S) the errors are typically larger by a factor 4. 

In Fig. 4 we plot kj vs. kq such that the masses determined for the Ising- 
and the Gaussian model at corresponding values of the reduced coupling 
T =\k,—k,c \ / i^c agree. On an infinite lattice one expects oc r/ | In rj- 1~3 and 
the extrapolation to the critical values kq^c thus should provide an estimate of 
c- On finite size however, the mass does not vanish at the phase transition. 
Still, within the accuracy of our results the points obey the scaling law quite 
accurately and we confirm that the leading critical exponents of m(r) agree 
in the two cases; the logarithmic factor cannot be resolved. 

Comparing the RI properties of our Ising model results for G{t,x) with 
those of the Gaussian version for same mass we find surprising agreement 



(cf. [0]). In Fig.s la and lb we compare the inverse Fourier transformed 
propagator G~^{p,ip) with the results due to the Gaussian model, to point 
out the good agreement. Indeed, within the error bounds we could not find 
any systematic deviations from the Gaussian behaviour. This conforms with 
the expectation that the boson propagator in the broken phase of the Ising 
model is in leading order that of a free particle like in the Gaussian model 
- up to logarithmic corrections. But for the rotational symmetry properties 
of the propagator and the considered quantities Af those seem to be small 
enough to permit neglecting them. This justifies to use the Gaussian values 
as our measure for violation of RI at given mass, at least for the mass values 
considered here. This is what we do in the subsequent discussion. 

In Table 2 we summarize our results for mass, / = ($_r) = {^)/ VZ, wave 
function renormalization constant determined from (126|), the ratio R of (0) 



and the value of N'g, one of our measures of violation of rotational symmetry 
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of the boson propagator. The large errors for the results S are due to the 
different updating method used there. In Fig.s 5 and 6 we plot rriR and / vs. 
the coupling and compare with the fit to the expected asymptotic behaviour 
m, ra for r ^ 



1 



mR ~ C^t 2 I In t 



/ ~ Cft^ |lnt|^ . (29) 



In Tab. 3 we give the results of our joint fit to the data for mn and /. The 
ratios C^/Cf are compatible with the asymptotic value 47r a/2/3 ||10| . 

Fig. 7 gives R vs. l/mn, where one could in principle read off the 
triviality bound. Comparing the results for actions C and S we find a 
difference similar to that observed in the recent study of 0(4) models in [|]. 
However, as discussed, it may be more appropriate to use another quantity 
than l/mpi, in order to compare the action (and regularization) dependence 
of that bound. We return to that point later. 

From these data one may obtain the constant C (cf. |jlO|) relating mn 
to the renormalized coupling, 

C = lim Imn {piXnf'/^'e^/^'^^] . (30) 

Depending on the closeness to the continuum limit, this fit to the asymptotic 
behaviour may be not very reliable. Also, C is determined just from our 
data for lattice size 16^ and should be extrapolated to infinite size in a finite 
size scaling study. Within these limitations we obtain the numbers given 



in Tab. 3. For action C Liischer and Weisz|]IO| determined the asymptotic 
value of C = 5.3(1.2); this is definitely larger than ours. However, fitting 
their results in the region of mass values considered here, we get a value of 
~ 3.0(3). We are obviously not asymptotic enough for C. 

As discussed it may not be wise to compare the values of R for different 
actions with each other at some lattice mass m^. Instead we follow Q 
and plot in Fig. 8 the renormalized coupling vs. A/'g, our measure for the 
amount of violation of RI. 

First we notice, that the difference between the results for C and C" and 
between those for F4 and S appears to become smaller, as compared to the 
presentation in Fig. 7; it is typically less than 10%. However, the difference 
between C, C on one hand and S on the other appears to be sizeable, 
about 30-40%. 
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In the results for the 0(4) model for F4 and S the authors es- 
timated the regularization dependence by comparing the results at given 
value of another measure of RI, derived from a perturbation expansion of the 
Goldstone-Goldstone scattering cross-section. There the difference between 
-F4, S compared to C was smaller than ours, of the order of 10-20 % for 
1-10% violation of RI. Although this variation could be due to the study of 
different models (one component vs. four component fields) this may also 
indicate a strong dependence on the choice of measure for RI. Maybe a sys- 
tematic study of higher order terms, like the one initiated recently [0 helps 
to clarify this issue. 

Let us summarize our results for the Ising model: 

• The similarities in the propagator between the Gaussian version and the 
Ising version within the same type of action (regularization) are bigger 
than between Gauss <-> Gauss or between Ising ^ Ising for different 
lattice actions. Therefore the RI behaviour for the discussed range of 
couplings seems to be dominated mainly by leading kinematical effects, 
which are already given in the corresponding Gaussian model. 

• The difference between the two measures of RI, Afc and A/'^ are less 
than 1%. This result is due to the close relationship between the cor- 
relation length and the momentum space propagator. 

• The regularization dependence of R comparing C with C and with 
S is less than 10% in the considered domain of couplings. 

• The regularization dependence of R comparing C or C with F4 or S 
is about 30-40 % in the considered domain of couplings. 

The last mentioned result allows two interpretations: The first is, that 
the regularization dependence of the envelope indeed may be about 40 %. A 
similar conclusion was put forward in ||2^ , |^ . On the other hand the studies 
[0, 1^ for the 0(4) model do not support this strong variation. So the second 
interpretation may be that there is substantial sensitivity on the definition 
of a measure of RI. 

Leaving aside the issue of the regularization dependence of the triviality 
bound, we found that different lattice actions lead to propagators with re- 
markable rotational symmetry, quite contrary to the naive expectation. The 
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leading violations are well described already by the corresponding Gaussian 
models. 
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Tables 

Table 1: The parameters for the four actions considered, defined by embed- 
ding into the hypercubic lattice. 



I 


N 




a 


{e,, i/ = l,...,iV} 


c 


8 


1 


1 




c 


32 


1 


7 


±e^, ±(e^ ± ex), X < n 




24 


1 


6 


±(e^±eA), \< IJ. 


s 


16 


771-^=1, ?79_16=-1/16 


3/4 





Table 2: The parameters for our results in the broken phase of the Ising 
limit. 



Action 








Zr 


R 


^fG 


C 


0.1507 


0.244(4) 


0.0832(1) 


0.929(2) 


2.93(4) 


7.6(2)E-4 




0.1513 


0.328(3) 


0.1043(6) 


0.924(4) 


2.14(4) 


1.37(2)E-3 




0.1520 


0.407(3) 


0.1227(5) 


0.918(6) 


2.31(3) 


2.09(3)E-3 




0.1525 


0.449(1) 


0.1339(2) 


0.909(3) 


2.35(1) 


2.55(l)E-3 




0.1550 


0.637(2) 


0.1769(1) 


0.886(2) 


2.60(1) 


5.07(3)E-3 




0.1560 


0.698(5) 


0.1914(3) 


0.873(3) 


2.65(3) 


6.06(8)E-3 


c 


0.0336 


0.231(1) 


0.0896(1) 


0.974(2) 


2.58(1) 


9.80(9)E-5 




0.0338 


0.319(2) 


0.1198(2) 


0.965(2) 


2.66(2) 


1.88(9)E-4 




0.0340 


0.393(6) 


0.1415(3) 


0.964(5) 


2.78(5) 


2.88(9)E-4 




0.0346 


0.560(10) 


0.1903(3) 


0.955(3) 


2.95(4) 


6.0(2)E-4 




0.0354 


0.738(4) 


0.2372(4) 


0.944(3) 


3.11(2) 


1.07(l)E-3 




0.0460 


0.255(3) 


0.0875(1) 


0.959(2) 


2.91(3) 


1.79(8)E-6 




0.0470 


0.474(1) 


0.1460(1) 


0.936(1) 


3.27(1) 


2.12(2)E-5 




0.0482 


0.670(2) 


0.1873(2) 


0.930(2) 


3.58(1) 


8.34(9)E-5 




0.0502 


0.917(4) 


0.2368(6) 


0.913(4) 


3.87(2) 


2.85(5)E-4 


s 


0.1653 


0.23(2) 


0.078(2) 


0.88(4) 


3.0(4) 


0.7(2)E-5 




0.1660 


0.33(1) 


0.100(2) 


0.89(2) 


3.3(2) 


3.0(4)E-5 




0.1675 


0.47(2) 


0.134(2) 


0.85(3) 


3.5(2) 


1.2(2)E-4 




0.1683 


0.55(3) 


0.147(4) 


0.85(4) 


3.7(3) 


2.3(5)E-4 




0.1694 


0.61(2) 


0.165(2) 


0.82(2) 


3.7(1) 


3.4(4)E-4 
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Table 3: Results for k^^^-*, C4 and Cf from a joint fit to m/j(r) and /(r) 
according to (pOj) ; we also give C according to a fit to (^). 



Action 


4^^^ 


C4 


Cf 


C 


C 


0.14984(2) 


4.24(1) 


0.636(1) 


2.9(1) 


c 


0.03340(1) 


3.76(1) 


0.674(1) 


6.4(2) 




0.04559(1) 


3.33(1) 


0.550(1) 


3.2(1) 


5 


0.16447(8) 


4.38(8) 


0.624(5) 


2.4(1) 
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Figures 



Fig. 1: (a) The inverse momentum space propagator G^^{p,(f) for the 
Gaussian model (action type C) for m = 0.407; the points give the numeric 
results from a simulation of the corresponding Ising model at coupling 0.1520 
such that m agrees (cf. Table 2). (b) The values of G~^{p = vr/2, ip) compared 
to the corresponding values of the Ising model simulation (resulting from an 
interpolation of points along the cartesian directions on the 16^ lattice grid). 

Fig. 2: Lines of constant values of C,~^{fn,ip), as derived in the Gaussian 
model for infinite lattice size, for the four actions considered. For small m 
(small lattice spacing) rotational symmetry is restored. 

Fig. 3: Our measures for the violation of RI A^g (full lines) and AT^ (dashed 
lines for m< 1. 

Fig. 4: We plot kj vs. such that the masses determined at corresponding 
values of r agree (action C; triangles: L = 8, squares: L = 16, circles: results 
for infinite size due to UTOl). 



Fig. 5: f vs. K, for the four actions considered; the curves give the fits to 
the expected critical behaviour (^). All data is for L = 16; in the figure for 



C the triangles denote the results for infinite lattice size due to [10 



Fig. 6: vs. k, for the four actions considered; the curves give the fits to 
the expected critical behaviour (p3). All data is for L = 16; in the figure for 



C the triangles denote the results for infinite lattice size due to [|10 



Fig. 7: R vs. l/rriR for the four actions considered (bars: C, squares: 
C, triangles: F4, asterisks: S); the curves give the fits to the asymptotic 
behaviour (IHDI). All data is for L = 16. 



Fig. 8: The ratio R vs. A^g; the curves denote fits to the asymptotic 
behaviour R{mR), where rriR is related to the corresponding value of Afc as 
discussed. The left-hand part gives the results on an extended scale. Symbols 
are like in Fig. 7. 
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